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(N ; ABSTRACT 

Context. Strings and other alternative theories describing the quantum properties of space-time suggest that space-time could present 

a foamy structure and also that, in certain cases, quantum gravity (QG) may manifest at energies much below the Planck scale. One 
\ of the observable effects could be the degradation of the diffraction images of distant sources. 

Aims. We searched for this degradation effect, caused by QG fluctuations, in the light of the farthest quasars (QSOs) observed by the 
■ Hubble Space Telescope with the aim of setting new limits on the fluctuations of the space-time foam and QG models. 

Methods. We developed a software that estimates and compares the phase variation in the interference patterns of the high-redshift 

1 1 1 QSOs, taken from the snapshot survey of HST-SDSS, with those of stars that are expected to not be affected by QG effects. We used 

a two-parameter function to determine, for each test star and QSO, the maximum of the diffraction pattern and to calculate the Strehl 
ratio. 

Results. Our results go far beyond those already present in the literature. By adopting the most conservative approach where the 
(2JT)| correction terms, that describe the possibility for space-time fluctuations cumulating across long distances and partially compensate 

for the effects of the phase variations, are taken into account. We exclude the random walk model and most of the holographic models 
of the space-time foam. Without considering these correction terms, all the main QG scenarios are excluded. Finally, our results show 
the absence of any directional dependence of QG effects and the validity of the cosmological principle with an independent method; 
that is, viewed on a large scale, the properties of the Universe are the same for all observers, including the effects of space-time 
, fluctuations. 

o . 

Key words. Gravitation, (Galaxies:) QSOs: general, Methods: data analysis, Methods: statistical, Cosmology: theory, Elementary 
\Q • particles. 
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y—i ' 1. Introduction ical phenomena and macroscopically observable effects, such 

7— I . as the production of mini-black holes at energies much lower 

1 Quantum gravity (QG) can be considered the Holy Grail of mod- than the classical Planck scale, like those expected in the large 

• i-h ! ern physics. Because of the general Heisenberg's uncertainty hadron collider (LHC, 2010; Dvali & Sibiryakov, 2008; Koch et 

^ • principle, in most of the approaches to the unification of grav- al., 2005; Hossenfelder, 2006; Hossenfelder, 2002). Other mod- 

^ [ itation and quantum physics, the classical geometrical concept els even describe the possibility of loosing of phase coherence 

of a continuous and smooth space-time support of general rela- in the light of distant sources, such as the farthest supernovae 

tivity is replaced, at very short distances, by a foamy geometric (SNe) and distant quasars (QSOs). Light, when propagating in 

scenario where fluctuations on the order of the Planck length the vacuum, will be affected by the fuzzy structure of space-time, 

l p =i 1.616 x 10~ 35 m occur at certain characteristic frequencies with the result of progressively blurring the diffraction images 

that are related to the Planck time t p a 5.391 x 10 44 s. For the of those sources (Amelino-Camelia, 1994; Amelino-Camelia, 

sake of simplicity, here we will not consider in this discussion 2000; Ng & van Dam, 2000; Ng & van Dam, 1994; Ng, 2003; 

more exotic formulations of Heisenberg's principle, whici can Christiansen et al., 2006). 
be found in the vast literature of quantum gravity. 

Without going into too much detail, each model of QG can be 

naively characterized by its peculiar space-time and energy flue- Naively, following the already cited literature, one expects 

tuation scale, which may occur also at frequencies that are dif- that a photon with energy E — hv is affected by an energy inde- 

ferent or even much lower than those expected at Planck scales, termination on the order of SE/E ~ dp/ p ~ (E/E p ) a , where a 

l/tp. The reader can find quite a vast literature describing all the is the characteristic exponent of the power-law distribution of a 

possible scenarios and behavior of gravity on quantum scales, phenomenological model that describes the energy fluctuations 

because up to now there has been no theory of QG. This scale induced in any quanta by QG effects because of the "grainy" 

is thought to be accompanied by the occurrence of new phys- structure of space-time. The quantity E p = h/tp is Planck's en- 
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ergy. The energy-momentum relationship at the lowest perturba- 
tion order is a function of the two parameters aq and a, 



p 2 = E\l±ao(E/E P yi 



(1) 



If this "grainy" structure of space-time emerges on time 
scales of f > tp, any time interval can be determined with an un- 
certainty t > f and intrinsic standard deviation derived from the 
following relationship, cr,/f = f(t p jf), where / > 1 for t < tp. 
When, instead, t » tp, one obtains / = (t p /t)(ao + a\{t p jt) + 
a 2(t p /t) 2 + ...) ~ ao(t p /t) <sc 1, where oq is the lowest order am- 
plitude of the quantum fluctuations in Planck units. This means 
that this foamy structure of space-time induces a random dis- 
persion of the photon arrival time of a bunch of photons emit- 
ted by a source. It is possible to reveal this effect by analyz- 
ing the phase coherence of the light of distant sources, as de- 
scribed by Lieu & Hillman (2003) and Ragazzoni et al. (2003) 



Equivalently, this dispersion is said to generate two differently 
fluctuating wave velocities that may systematically cumulate 
during the journey of each photon. The wavelength A of the trav- 
eling photon naturally defines the minimum and characteristic 
scale length on which physical quantities such as, the phase ve- 
locity (v ; ,) and the group velocity (v f ), and their respective dis- 
persion relationships can be defined. 

The phase of the electromagnetic wave, associated with the 
photon propagation, after having crossed the distance A in the 
tim e t = A/v,,, beco mes </> = 2-nv p jv R . |Lieu & Hillman (2 003) 



and Ng et al. (200 1 ) claim that QG fluctuations may induce well- 
defined phase indeterminations such that an electromagnetic 
wave, traveling a long distance L, may cumulate a phase vari- 
ation on the order of 



(A(/>Qc)n 



l a ]_} a 

2nat)— 



2n- 



1 ± V2a 



he 
AE~ P 



(2) 



if the emission and detection wavelengths remain equal during 
the photon propagation. The parameter c is the speed of light and 
ao ~ I (Ng, 2003). If, instead, one also considers the effect of 
cosmological redshift z, a larger phase dispersion is obtained, as 
proposed by Steinbring (2007) 



(A(p QC ) z ~ 2na Q 

x f(l+z)L- a 
Jo 



l a P {\ -a)c 
A Hoqo 

l - qo 



1 - 



dz, 



(3) 



where Hq is the Hubble constant, qo the deceleration parameter 
of the Universe, and the luminosity distance is given by 



L = 



[q Z - (1 - qo)( V 1 +2<?oZ - 1)] 



(4) 



This oversimplified phenomenological description of QG 
models the foamy structure of space-time with only the two pa- 
rameters flo and a at the first order in amplitude fluctuations. 
By varying a, one obtains the three main scenarios of QG. The 
random walk model has a = 1/2 and random perturbations of 
space-time add incoherently with a square root dependence in 
the simplest case, or even following a more general distribution 
function. Examples are the works by |Amelino-Camelia~2 000) 
and Ng & van Dam (1994) The holographic model, which is 
consistent with Beckenstein-Hawking's holographic principle 
and black hole entropy, has a = 2/3 ('t Hooft, 2003). The value 



a — 1 describes the standard model of QG in which the gravi- 
tational field fluctuations are close to the Planck time (Smolin, 
2003; Rovelli, 2004). 

Those effects can in principle be detected by analyzing pos- 
sible modifications of the interference fringe structures using 
stellar interferometry. Equivalently, the cumulated phase vari- 
ation can also be detected by imaging and analyzing, with a 
diffraction-limited telescope, the ring structure of Airy disks of 
cosmological point-like sources (Ragazzoni et al., 2003). Airy 
rings are still observed if the cumulated phase is limited to the 
value 6(p < 2n. 

Previous results, based on data from stars and extragalac- 
tic sources (Ng & van Dam, 1994; Lieu & Hillman, 2003), 
seemed to rule out both the random walk model and the holo- 
graphic principle and to put severe limits on the standard model. 
Similarly, the standard model also appeared to be ruled out by 
data on PKS 141 3+35 (z=0.247) and by the detection of the Airy 
rings from SN 1994D (z= 0.001494) (Ragazzoni et al., 2003). 
However, in these conclusions the authors did not consider the 
more recent (Ng, 2003) possibility that the space-time fluctua- 
tions may cumulate across long distances and partially compen- 
sate for the effects of the phase variations. In this case the correc- 
tion factor is (L/A)~ a and only the random walk model could be 
excluded with the data available so far. A more sensitive anal- 
ysis performed by Steinbring (2007), about the Strehl ratios of 
high-redshift QSOs, pointed out a possible deterioration of the 
interference structures due to the space time foam. 

In this paper we resume the test on QSOs with a new ap- 
proach to the data analysis and we put new upper limits on the 
QG fluctuations by analyzing the images of distant QSOs, ob- 
served by the Hubble Space Telescope (HST) with a software 
ad-hoc developed for this purpose. 

With our estimates, if we neglect those possible cumulat- 
ing phase effects, all the models of space-time foam would be 
immediately ruled out, thereby posing strict limits on the stan- 
dard model. However, in our paper we use a more conservative 
approach, in which the space-time fluctuations are expected to 
compensate for their variations owing to their "averaging" and 
"smoothing" during the propagation of the light beam, causing 
the lack of detection of space-time fluctuations. 



2. Data analysis 

We searched for QG effects on the interference patterns of the 
high-redshift QSOs taken from the archival HST ACS snapshot 
survey of Sloan Digital Sky Survey (SDSS) sources (Proposal: 
9472, PI: Strauss) and described by |Richards et al. (2006)1 The 
dataset consists of 157 sources located in the redshift interval 
4 < z < 5.4 and acquired with filter F775W (A c = 776 .4 
nm, bandwidth AA = 152.8 nm). Four other sources, located 
at redshifts 5.7 < z < 6.3, were instead taken with the filter 
F850LP (A c = 944.5 nm, bandwidth = 122.9 nm). The de- 
tector used in the observations was the High Resolution Camera 
(HRC), which has a pixel scale of 0.0246 arcseconds. We used in 
our analysis the drizzled (DRZ) images that were corrected for 
geometric effects. The DRZ files were all reduced by adopting 
the standard procedure, including the overscan-, bias-, and dark- 
corrections and the flat-field and the photometrical calibration. 

The luminosity distances of the QSOs, L, were calcu- 
lated by assuming a concordance Universe model described in 
|Krauss (2002)1 H = 72 km/s/Mpc, (fi M , C2a)=(0.3, 0.7), and 
also by taking the small effects on A due to the expansion of the 
Universe into account. 
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To determine the QG effects through the progressive image 
degradation of a distant source, we estimated the Strehl ratio S of 
the source for a given wavelength A through the empirical phase 
variation 



SDSS J015339. 60-001104. 8 - F775W 



A 



- V-MS) 



(5) 



where D = 2.4 m is the telescope diameter. The Strehl ratio S is 
defined as the ratio between the actual observed peak intensity 
(Max) source and what is theoretically expected from an unaber- 
rated telescope (Max)th 



(Max) s , 



(Max),, 



exp 



(6) 



From the Strehl ratio S we obtain the total phase variation A<p 
that for QSOs can be written as 



A<£(z)qso = A0 ab + A<p(z) siz e + A0(z) QG + A0(z) kns . 



(7) 



where A0 a b is the phase variation due to the aberrations present 
in the optical path of HST, A0qg those due to QG, A0 s ; ze to the 
resolvability of QSOs and/or their host galaxies, and A0(z)i ens to 
the gravitational lensing. The phase variation due to gravitational 
lensing, averaged in all the exposures of the dataset, is expected 
to behave as an increase in the blurring already caused by the 
QG effects, so it cannot, in principle, be easily separated from 
the latter ones except in special cases. The presence of spikes, 
caustics, and other typical effects of microlensing were not been 
observed in our sample. Because of this, the detection of QG 
effects would be overestimated, because unavoidably entangled 
with gravitational lensing and optical aberration effects. 

To reduce the blurring caused by the aberrations of the tele- 
scope, we modeled this image degradation effect by analyzing a 
sample of pictures of stars. The effects of QG are negligible at 
non cosmological distances and A0 stal - A0 a b. 

In the case of the QSOs, we now consider the quantity 



(Acf>(z) Q so) Corrected 



A0(z) QSO - <A0) star 

A0(z) size + A0(z) QG + A0(z)i ens : 



To adopt a more conservative approach, we use the averaged 
value of the phase variation (A<f>) stai instead of the dispersion 
of A0 stal values found in the sample of the stars analyzed by 
us to calibrate our procedure. The quantity 6<f) mrs includes the 
errors due to the data analysis, the different positions on the de- 
tector and the different colors of the stars. The indeterminant 
64- 'color ^ 10~ 8 rad represents the error due to the different colors 
found between the stars and QSOs. This quantity has been es- 
timated by analyzing the PSFs generated by the Tiny Tim soft- 
ware. 

Our approach is different from the procedure applied by 
Steinbring (2007) in which the Strehl ratios of QSOs were not 
corrected but only compared to an appropriate PSF generated for 
each of the QSOs by the Tiny Tim software. In fact, we did not 
assume any particular PSF model in our investigation, and we 
corrected the phase variation of QSOs by using the correction 
terms derived empirically from galactic stars in Eq.©. 

For our analysis we developed a dedicated software to deter- 
mine the position and the maximum of QSOs and stars. For each 
source imaged by the ACS camera, we assumed that the central 
region of the interference structure has an elliptical symmetry. 
The position and the value of the maximum of the light pattern 
are determined in a region of 3 X 3 pixels. The light intensity 
distribution is then interpolated by a 2-D function 





HST IMAGE 



RECONSTRUCTED IMAGE 



RESIDUALS 



Fig. 1. Left: HST image of the QSO QSO 0153-001 1, IAU name 
SDSS J015339.60-001 104.8 located at z = 4.205. The field of 
view is 15 x 15 pixels. Center: we show the central region (3x3 
pixels) reconstructed using the function F (EqO- We plot the 
cartesian axes centered on the maximum of the intensity pattern. 
Here e = 0.23 and a — 0.14 rad. Each pixel of the analyzed area 
is divided into 15 x 15 sub-pixels. Right: we show the residuals 
of the central analyzed area (3x3 pixels). 



F (x,y, [mx, cx, cy, e, 8, cr\, era]) = 



mx 



x ) ^ <x 3 _„ exp 

n=\ 



x\(x, • • •) - ex\(x, ■ ■ •) + y\(x, ■ ■ ■) 



2^(1 - e) 



(9) 



where the pair (cx, cy) identifies the position of the maximum 
(mx) in the Cartesian coordinate system, the origin of which is 
placed on the center of the brighter pixel. The quantity e is the 
eccentricity of the ellipse, and 9 is the rotation angle between 
the axes of the ellipse and the cartesian axes. The parameters cr\ 
and o~2 modify the shape of the function F, and we write the 
argument of Eq.|9]as 



x\(x,y, cx, cy, 0) 
y\(x,y, cx, cy, 9) 



(x - cx) cos 0—(y— cy) sinfl 
(x - cx) sin# + (y - cy) cos 9 . 



(10) 



With a recursive procedure, we tune the values of the param- 
eters mx, cx, cy, e, 9, cr, , cr 2 of the 2D function, F, until minimiz- 
ing the difference between the measured and the fitted intensity 
of each pixel, on a 3 x 3 matrix situated around the brightest 
pixel. In the calculations, each pixel is divided into 15 x 15 sub- 
pixels, and for each of these divisions the procedure determines 
the intensity as a function of the position in the detector. At the 
end, we obtained both the position of the center and the intensity 
of the light pattern in the 3x3 area with a precision of 15 x 15 
division for each pixel (see Figs[T]and|2]i. We found deviations 
between the fitted and measured intensities to be less than 1%. It 
is important to point out that we considered only the light distri- 
bution inside the region delimited by the inner Airy diffraction 
ring. The values of the maxima of stars and QSOs are normalized 
with respect to the flux inside an area, centered on the brightest 
pixel, with a radius of 2 arcseconds. An example is reported in 
Figs. 1 and 2 with the QSO QSO SDSS J015339.60-001 104.8 
(IAU) (z=4.205). 

In Fig. [3] we plot the phase variation of a sample of stars 
chosen with the values of the maximum of the diffraction pattern 
in the same range as the QSOs. We indicate with (A<p star ) the 
averaged value of the phase variation of the sample at which 
we associate the error 6<p srar defined as the maximum difference 
between the phase variation values of the stars and (A(f> star ). 

In Fig.|4]we plot the corrected phase variation of the QSOs, 
(A0(z)Q5o)cortBoted> at which we associate the error Stp st ar- We 
found no blurring effects of QGs. The phase variations observed 
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Phase Variation vs Redshift of Point-like Quasars' classes 
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- Measured Intensity 

- Intensity of the reconstructed image 
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Measured Intensity 

Intensity of the reconstructed image 
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Fig. 2. Top: Id cut along the x-axis (solid line) and the y— 
axis (dashed line) of the measured intensity of QSO SDSS 
J015339.60-001 104.8. We analyzed the area within the two dot- 
ted lines. Bottom: we show the measured intensity (red line), the 
interpolated function F (solid blue line), the intensity of the re- 
constructed image (dashed blu line), and the residuals (dotted 
blue line) of the 1-d cut along the x-axis (left) and the y-axis 
(right) of the central area (3x3 pixels). 
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Fig. 3. Phase variation of 72 stars imaged with the High 
Resolution Camera with the filter F775W (upper panel) and 36 
stars with the filter F850LP (bottom panel). We plot also the av- 
eraged value (A(p star ) (solid line) and the error 6tp star . The dashed 
lines show the values (A0 J/ar ) ± 6<f> star - 



are much smaller than those expected from the space-time fluc- 
tuations described by the holographic principle model having 
a - 2/3 and a Q = 1. The phase variation (A0(z) eso ) Co iTected 
and the redshift z, follow, at the first order, a linear relation- 
ship (A0(z)Qso)conected = az + b, whose coefficients are a = 
(0.26 + 0.78) x 10~ 8 and b = (-0.5 ± 3.5) x 10~ 8 , indicating 
the possible resolvability of QSOs and/or their host galaxies. We 
make use of the relations 



\— and a » 



(ID 



to evaluate the errors on the coefficients a and b, where N is the 
total number of QSOs and A = N £ z- - (£ z*) 2 - 

To put our limits to the two parameters of the QG, a and 
a , we assume that (A0(z) eso )c O nected > A(p(z) QC - Our limits 
are conservative because in the quantity {A^>(z)qs o)corrected we 




Fig. 4. Corrected phase variations of QSOs as a function of the 
redshift z. No blurring effects due to QG are observed. The lin- 
ear interpolation of all the data sets show a positive trend, indi- 
cating the possible resolvability of the sources - blue solid line 
(color online). The red (color online) continuous and dashed "z"- 
shaped lines represent the expected growth of the phase varia- 
tion, A(p, as a function of the redshift z- The break of the two 
curves indicating the limits of QG effects with (red solid line) 
and without cosmological redshift corrections (red dashed line), 
are due to the change of filters, F775W and F850LP. 



cannot disentangle the QG effects from those from the resolv- 
ability of the QSOs. In Fig. [5] we compare our results with the 
upper estimates found in the literature, which refer to the esti- 
mates obtained with sources located at distances up to z — 5.34 
(Ragazzoni et al., 2003). 

As reported in Fig. [6] our results exclude the random walk 
model and most of the holographic scenario of the space-time 
foam. For a - 1, one obtains (A0(z) GlSO ) Co iTected = (0.26 z - 
0.5) x 10~ 8 rad. By using the results previously discussed and the 
errors on the estimate of the coefficients a and b, one obtains the 
lower limit a„,j„ = 0.677, with an 83.3% CL lower limit without 
considering the expansion of the universe, and a z = 0.684 with 
an 85.5% CL lower limit, if we take the cosmological redshift of 
the photons emitted by QSOs to the observer's wavelength into 
account. Both the results exclude the Holographic scenario, but 
for low values of the the lowest-order amplitude of the quantum 
fluctuations, ao < 5x 10~ 2 and ao < 0.15. The expected value for 
the amplitude of the quantum fluctuation is ao ~ 1 (Ng, 2003). 

3. Conclusions 

From the analysis of the diffraction images of the most distant 
QSOs, obtained with HST, we found no blurring effect that could 
be caused by the interaction of photons with quantum grav- 
ity fluctuations. Our results indicate that most of the models of 
space-time foam consistent with the holographic principle and, 
consequently, with black hole entropy, can be ruled out for most 
of the values of the fluctuation amplitude parameter ao- 

Surprisingly, we did not find any indirect trace of the mag- 
nification of QSO images expected from the geometry of the 
Universe adopted in our calculations that could generate a possi- 
ble additional blurring in the diffraction images combined with 
that expected from QG effects. The linear regression found in 
our estimates indicates that an eventual combination of QG ef- 
fects and that possible of QSO magnification does not show any 
detectable contribution with the dataset available. If a significa- 
tive trend were found, we should have estimated the QSO lumi- 
nosities and calculated a QSO luminosity function (Richards et 
al., 2006) to disentangle this contribution from the one of QG 
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Fig. 5. Parameter space ao and a The red region (colour online) 
in the left is the portion of the parameter space that is discarded 
by our observations with the most conservative approach that 
takes in consideration the cumulating effects of space-time fluc- 
tuations (see text). The white region in the middle is the addi- 
tional portion of the parameter space that would be discarded 
if no compensation due to the cumulating effect of space-time 
occurs. The grey region in the right is the portion of the pa- 
rameter space that cannot be excluded in any case. The data 
obtained with our analysis from the farthest QSOs exclude the 
random walk and the holographic models of space-time foam in 
the range of ao, here reported. These results, indicated by the la- 
bel QSOs, go far beyond the previous estimates reported in the 
literature (indicated by z=5.34 Universe), obtained without con- 
sidering the cumulating effects of space-time. 













Holographic Model 
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Random Walk Model 
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Fig. 6. Limitations to the parameters ao and a from QSOs. The 
excluded regions are below the dotted-dashed or the solid line 
depending on whether the photons redshifted to the observer 
wavelenghts, were taken into account. The holographic model 
(a = 2/3) of space-time foam is preserved for ao < 5 1CT 2 and 
ao 55 0.15 (upper panel). The random walk model (a = 1/2) is 
preserved for ao 55 10~ 13 and ao 55 10~ 12 (lower panel). 



fluctuations. One can suppose that either the Universe is flat in 
the interval 4 < z < 6 or that QSOs may present some unknown 
evolutionary effects at large z. Future observations with EELT, 
when equipped with adaptive optics, of QSOs at higher redshift 
also imaged with a filter at lower wavelengths should improve 
our results. 



Finally, our results show the absence of any directional de- 
pendence on QG effects, since QSOs are located in very different 
regions of the sky at redshifts between 4 < z < 6.3. This would 
also limit the effects of intergalactic medium on the light across 
these cosmological distances in those observation bands. This 
confirms the cosmological principle on large scales. The isotropy 
on large scales of the Universe, verified by the distribution of su- 
perclusters of galaxies and QSOs on very large scales, is also 
verified through the independence of the behavior on smaller 
scales of the quantum fluctuations from any direction in space. 
Viewed on a large scale, the properties of the Universe and the 
laws of physics are the same for all observers. 

3.1. The effect of QG fluctuations in quantum experiments 

We can also infer from our results that quantum states will be 
preserved for very long distances because they are not affected 
by the blurring of space-time fluctuations. This is a very impor- 
tant result for quantum information and the foundations of quan- 
tum physics: a photon can travel undisturbed for billions of light 
years within the errors dictated by current HST technology. One 
can refer to Kok & Yurtsever (2003) Bouwmeest er et al. (2 000) 
and [Bertlmann et al. (2002) for deeper insight. With our results 
we can set an upper limit to the phase factor imposed by the 
foamy structure of space-time to the photons that travel from the 
farthest QSOs caused by QG fluctuations, which is on the order 
of50~2 1O- 8 . 

We consider an ideal quantum communication experiment, 
where Alice on the Earth wants to exchange single quanta of 
light with Bob located on the most distant QSO observed by 
HST, with wavelength A — 1 fim and build a single-photon cryp- 
tographic quantum key. The two parties will find that the proba- 
bility of having one qubit mismatch caused by quantum gravity 
is almost negligible with respect to the errors usually present in 
those experiments, affecting at maximum only one qubit over a 
set of ~ 5 x 10 7 qubits used in the building of the key. Earth- 
Moon and Earth-Mars experiments would cumulate a phase 
shift A(f> ~ 4.18 x 10~ 14 and A</> ~ 2.15 x 10~ 13 , respectively. The 
value a — 1 suggests a lower limit A(p ~ 6.15 x 10~ 17 when non- 
cumulating QG fluctuations occur on the Planck scales. Similar 
argumentations are valid for entangled photons. 
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